arXiv:1509.04552v3 [math.PR] 13 Jul 2016 


FIXED POINTS AND CYCLE STRUCTURE OF RANDOM PERMUTATIONS 


SUMIT MUKHERJEE 


Abstract. Using the recently developed notion of permutation limits this paper derives the lim¬ 
iting distribution of the number of fixed points and cycle structure for any convergent sequence of 
random permutations, under mild regularity conditions. In particular this covers random permu¬ 
tations generated from Mallows Model with Kendall’s Tau, ^ random permutations introduced in 
[11], as well as a class of exponential families introduced in [15]. 


1. Introduction 

Study of random permutations is an area of classical interest in the intersection of Combinatorics 
and Probability theory. Permutation statistics of interest is indeed a long list which includes number 
of hxed points, cycle structure, length of longest increasing sub-sequence, number of descents, 
number of cycles, number of inversions, order of a permutation, etc. Most of this literature focuses 
on the case where the permutation tt^ is chosen uniformly at random from Sn- For example it is 
well known that the number of fixed points of a uniformly random permutation converges to Poi{l) 
in distribution. More generally, denoting the number of cycles of length I by Cn{l), we have 

{^(1), • • • , Cnil)} 4 {Poiil), Poi{l/2), • • • , Poiil/l)}, 

where the limiting Poisson variables are mutually independent. However, not much is known in this 
regard outside the realm of the uniform measure. Possibly the most widely studied non uniform 
probability measure on Sn is the Mallows model with Kendall’s Tau, first introduced by Mallows 
in [13], which has a p.m.f. of the form 

Mn,,{7r) = (1.1) 

^n,q 

Here Inv{'Kn) '■= < 0} is the number of inversions in vr^, and 

g > 0 is a scalar parameter. In this case an exact formula is known for the normalizing constant 
Zn^q, and expectation and variance formulas for Inv{'Kn) are easy to derive (see for e.g. [9]). In [6] 
Borodin et al. asked the question of behavior of permutation statistics such as cycle structure and 
longest increasing sub-sequence for general class of Mallows models which includes the Mallows 
model with Kendall’s Tau. This question was partially answered by Mueller-Starr in [14], where 
they derived the weak law of the length of the longest increasing sub-sequence. Specihcally, for the 
scaling ra(l — q{n)) —>■ f3 they showed that 

^L/5(7r„) 4 £(/3), 

vn 
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where C{j3) := 2/3 ^/^sinh ^(Ve^ — 1) for /3 > 0, and 2|/3| ^/^sin ^(Vl — e^) for /3 < 0. For the 
scaling when n(l — q{n)) oo, it was shown by Bhatnagar-Peled ([5]) that 


1 

n^l — q{n) 


LIS{7Tn) 4 1 . 


The more recent work of Basu-Bhatnagar ([2]) consider the case q{n) = g 7 ^ 1 is fixed, and prove 
a weak law for LIS{'nn) (they also derive a central limit theorem for (/<!). This answers the 
question of LIS for Mallows model with Kendall’s Tau for all parameter scalings, at least at the 
level of weak limits. On the other hand, the question of the cycle structure still remains largely 
unanswered. See however the recent work of Gladkich-Peled, who derive the order of expected 
number of cycles in a Mallows random permutation in [10, Theorem 1.1], when the underlying 
parameter q{n) G (0,1) is arbitrary. 

In a different direction, in [11] the authors Hoppen et al. proposed a framework where a permu¬ 
tation can be viewed as a measure. This is described below in brief: 

For a permutation 7 r„ G Sn define the measure on [0,1]^ as 


1 

n 


n 

%/n,7r„(j)/n)- 

2=1 


A sequence of permutations {7rn}n>i with Tin G Sn is said to converge to a measure /x, if the sequence 
of probability measures r'jr„ converge weakly to //. Any such limit is in AI, the set of probability 
distribution on the unit square with uniform marginals. Any /r G Ad is called a permuton (following 
[11]), and it is shown in [11, Theorem 1.6] that any /x G AI can indeed arise as a limit of a sequence 
of permutations in this manner. See [3, 11] for a more detailed introduction to permutation limits. 

If {vr„}n>i is a sequence of random permutations (not necessarily in the same probability space), 
the sequence is said to converge to a deterministic measure /x G Al in probability, if the sequence 
of measures converge weakly to the measure /x in probability. Equivalently, for any continuous 
function / on the unit square, one has 


lim 

n^oo n \n 
2=1 


7rn(i) \ 
n J 



f{x,y)dn. 


Using the topology of permutation limits in [15] the author gave a new proof for a large deviation 
principle (originally proved in [19]), and used it to analyze a class of exponential families on the 
space of permutations. The large deviation principle was re-derived in [12], where Kenyon et 
al. study permutation ensembles constrained to have fixed densities of finite number of patterns. 

It was shown by Starr in [17] that if is generated from a Mallows model with Kendall’s Tau 
with parameter q{n) such that n(l — q{n)) —)> /3, then the sequence of measures converge weakly 
in probability to a measure G Al induced by the density 


. ^ __ (/3/2)sinh(/3/2) __ 

eA4 cosh(/3(x — y)/2) — e“A4 cosh(/3(x-|-y — l)/2) ’ 


( 1 . 2 ) 


which is the Frank’s Copula (see [16]). Since 7 r„ converges weakly to the measure /Xp^, in an 
attempt to understand the marginal distribution of 7 r„(x) one might conjecture that P„( 7 r„(x) = 
j) ^pp{i/n,j/n). We will show that this is indeed true, under certain regularity of the law of 
the random permutations. We start by introducing some notations. 


Definition 1.1. For I G [n] := {1, 2, • • • , n} let 

S{n,l) := {p := {pi,P 2 , ■ ■ ■ ,Pi) G [nf : p^ ^ pb for all a^b,a,b£ [/]}. 
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Then we have \S{n,l)\ = For p,q G S{n,l) let ||p — q||oo := \pa — qa\- Also for 

p G 5(n, 1) let vr„(p) denote the vector (7rn(pi), • • • , TTn{Pk))- 


For every n > 1 let iTn be a random permutation on Sn with law In [3, Def 6.2] the authors 
define a notion of equi-continuity of random permutations, which they show is implied by the 
condition 


lim lim sup 

(5—>0 n—>-oo p^q.rg5(n,i):||p-r||cKj<n(5 


IPn(7rn(p) = q) 
IPn(7rn(r) = q) 


(1.3) 


(see [3, Prop 6.2]). In particular for I = 1 condition (1.3) in spirit demands that the function 
IF’n('Tn(p) = q) is equi-continuous in p. In this paper we will need an extra notion of equi-continuity 
which demands that the function P„(7r„(p) = q) is jointly equi-continuous in p,q. This is stated 
below: 


Definition 1.2. A sequence of random permutations TTn is said to be equi-continuous in both 
co-ordinates if 


lim lim sup 

5—>-0 n—>-oo p^q_r^secS(n,Z):||p—r||cx><n<5,||q-s||<n(5 


IP’n(7rn(p) 

IF’n(vrn(r) 



= 0 . 


(1.4) 


Definition 1.3. Let C denote the set of all strictly positive continuous functions p on [0,1]^ with 
uniform marginals, i.e. 

/ p{x,y)dx= p{x,y)dy = l. 

Jo Jo 

Denote by //p G Af the measure induced by p. 

Our first theorem now proves an estimate of Pn,(7r„(p) = q) for vectors p, q if is equi- 
continuous in both co-ordinates, and converges in the sense of permutation limits to pp. 


Theorem 1.1. Suppose {'Kn}n>i is a sequence of random permutations with Tin G Sn, such that 
the sequence is equi-continuous in both co-ordinates, i.e. it satisfies (1.4). If{TTn}n>i converges to 
Pp for some p G C, we have 


lim sup 

p,qe5(n,Z) 


n^Pn(7rn(p) = q) 

nL, (>(??. s) 


(1.5) 


As an immediate corollary of Theorem 1.1 we obtain limiting distribution of the vector 7rn(p). 
A more general version of this corollary was already derived in [3, Proposition 6.1]. 


Corollary 1.2. Suppose pn G S{n,l) is such that 

lim -p„ = X G [0,1]^ 

n^oo 71 

//{vTn} n>i is a sequence of random permutations with Tin G Sn which satisfies (1.5) for some p G C, 
then 

-7rn(Pn) {i^(a^i), • • • 
n 

where {T(xa)}a=i '^xe mutually independent with Y{xa) having the density p{xa ,.). 


Having proved Theorem 1.1 we now turn our focus on the number of fixed points, or more 
generally the statistic 

n 

Nn{TTn, CTn) ■= ^ l{Tn(i) = CTnii)} 
i=l 
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for any cTn G Sn, which denotes the number of overlaps between 7r„ and ct„. In this notation the 
number of fixed points of equals A^('/r„, e„), where e„ is the identity permutation in Sn- By (1.5) 
^niT^riiC^n) IS approximately the sum of n independent variables, and so should be approximately 
distributed as Poisson. Our next theorem confirms this conjecture, showing convergence to Poisson 
distribution of A^„(7r„,(T„) in distribution and in moments. 

Theorem 1.3. Suppose {'Kn}n>i is a sequence of random permutations with ttu £ Sn which satisfies 
(1.5) for some p ^ C. If an converges to p,, then limn^oo'^IIn{T^n,an)^ = ^Poi{p[p])’^ for any 
A: G N, where p[p\ := ^2 p{x,y)dp, and Poi{\) is the Poisson distribution with parameter A. In 

particular this implies Nn{TTn,an) Poi{p[p]). 

Remark 1.1. Setting an = it follows by Theorem 1.3 that the number of hxed points in has 
a limiting Poisson distribution with mean p(x, x)dx, provided {vr„}n>i satishes (1.5) for some 
p£C. 

The random variable Nni^n, Sn) = XlILi I{'^n{i) = *} is essentially the number of cycles of length 
1, and a similar intuition for Poisson approximation holds for cycles of length I for any / > 1. In 
order to make this precise, we introduce a few more notations. 

Definition 1.4. For any I G [n] setting U{n,l) := {p G S{n,l) : pi = min(p(j,a G [1])} note that 
U{n,l) C S{n,l), and \S{n,l)\ = I x \U{n,l)\. For p G S{n,l) let r(p) G S{n,l) denote the vector 
{P2,P3, ■ ■ ■ ,PhPl)- 

As an example if / = 3 and n = 6 then the vector p = (2, 5,4) G U{n, 1), as 2 = min(2, 5,4). In 
this case r(p) = (5,4,2) G S{n,l) but does not belong to U{n,l), as 5 min(2,5,4). Thus T is 

the shift operator which shifts every co-ordinate by 1. 

For any / > 1 let 

Cn{l) ■■= = ^(p)} = 7 li^^(p) = ^(p)}- 

pGU{n,l) pGS{n,l) 

Then Cn{l) is the number of cycles of length I, where the factor I in the second definition accounts 
for the fact that every cycle is counted I times in the second sum. In particular we have C'n(l) = 
Uni'^rii e-n) to be the number of hxed points. Also let 

Cp(0 ■=t[ p{xi,X 2 ) • • • , p{xi,xi)dxi ■ ■ ■ dxi- 

The following theorem derives the limiting distribution for Cnil) under condition (1.5). 

Theorem 1.4. Suppose {TTn}n>i is a sequence of random permutations with iTn G S/n which 
satisfies (1.5) for some p £ C. Then for any {ki, - ■ ■ , A:;} G N* we have 

i i 

lim E TT Cnia)’^^ = TT EPoi(cp(o))^“. 

a=l a=l 

In particular this implies 

{^(1), • • • , Cn{l)} 4 [Poi{Cp{l)), • • • , Poi{Cp{l))], 
where {Poi{cp{i))'\\^i are mutually independent. 
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Remark 1.2. Thus the number of cycles of length I has a limiting Poisson distribution with 
parameter Cp{l), whenever the sequence of permutations satishes (1.5) for some p £ C. In 
particular if 7r„ is uniformly random then (1.5) holds for the function p = 1, in which case Cp{l) = j 
for alH > 1. In this case we get back the classical result that the number of cycles of length I is 
asymptotically Poi{l/l), and the random variables {(^^(l), • • • , Cn{l)} are mutually asymptotically 
independent for any I £ N. 

1.1. Applications. As applications of Theorem 1.3 and Theorem 1.4, we will now derive the limit 
distributions of the number of fixed points and cycle structures for three classes of non uniform 
distributions on Sn- 

(i) The hrst result in this direction is the next corollary, which deals with the Mallows model 
with Kendall’s Tau. 

Corollary 1.5. Suppose 7r„ is a random permutation on Sn generated from the Mallows 
model with Kendall’s Tau defined in (1.1), such that n(l — q{n)) -£■ fi £ (— 00 , 00 ). In this 
case the following conclusions hold with pp as defined in (1.2). 

(a) If {(Tn}n>i is a sequence of non random permutations with Un £ Sn converging to p, 
then Nn{Trn,crn) converges to Poi(/u[pp]) in distribution and in moments. 

(b) |c'„(l),--- ,Cnil)} converges to |Poi(cp^ (1)), • • • , Poi{cpp{l))'^ in distribution and in 
moments, where {Poi{cpp{i))}\^.^ are mutually independent. 

As an illustration of the Poisson approximation, in hgure 1 we compare the histogram 
of the number of hxed points in a permutation of size n = 100 with the limiting Poisson 
prediction. We used 10000 independent observations from the Mallows model with Kendall’s 
Tau with parameter q{n) = From the picture it seems that the Poisson prediction 

is fairly accurate for n = 100. Since q{n) < 1 it is expected that this model will have more 
hxed points than a uniformly random permutation, which is rehected in the fact that the 
mean of the Poisson distribution is much larger than 1. 

(ii) Another class of non uniform measures on permutations introduced by the author in [15] is 
the following: 

For any continuous function / on the unit square, let be a one parameter exponential 
family with sufficient statistic 

i=l 

More precisely, the p.m.f. is given by 

= ( 1 . 6 ) 

where Zn{f, 0) is the log normalizing constant of the model. In particular the permutation 
model obtained by the following two specihc choices have been studied in the Statistics 
literature: 

(a) f{x,y) = \x — y\, which gives the statistic ^11=1 K(^) “ ^1 known as the Spearman’s 
Footrule. 

(b) f{x,y) = (x —y)^, which gives the statistic X]r=i(^(*) known as Spearman’s rank 
correlation Statistic. 

See [8, Chapter 5,6] for more on these and other non uniform permutation models con¬ 
sidered in the Statistics literature. The convergence of a sequence of random permutations 
Tin generated from Qn,f,e of (1.6) was shown in [15, Theorem 1.4]. Building on this result. 
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Figure 1. Bar plot of empirical distribution (from 10000 observations) of number of fixed 
points in a permutation of size n = 100 from the Mallows model with Kendall’s Tau with 
parameter in green, compared to the Poisson prediction in yellow. 

the next corollary derives the limiting distributions of the number of fixed points and cycle 
structure for a permutation vr^ generated from this model. 

Corollary 1.6. Suppose -Kn is a random permutation on Sn generated from the model Qn,f,e 
defined in (1.6) for some function f which is continuous on the unit square. In this case 
the following conclusions hold: 

(a) The sequence {7r„}„>i converges weakly to a non random measure £ Ai with a 
continuous density gj 

(h) If {an} n>i is a sequence of non random permutations with Un £ Sn converging to p, 
then Nn{Trn,an) converges to Poi{p[gfp]) in distribution and in moments. 

(c) |c'„(l),--- ,Cnil)} converges to ^Poi{cg^ g(l)), ■ ■ ■ , Poi{cg^ g{l))'^ in distribution and 
in moments, where {Poi{cg^ g{i))}\^i are mutually independent. 

(hi) The final class of permutation models that we consider is a non parametric model with a 
measure as the parameter, as opposed to the previous two models which are one parameter 
models. This class of models will be referred to as p random permutations, and was hrst 
introduced in [11]. 

Given any p £ Ai let (Xi, Y)), • • • , {Xn,Yn) be i.i.d. random vectors with law p. Dehne 
a permutation G Sn as follows: 

If there exists I £ [n] such that Xi = X(^^pYi = Y(^jp then set TTn{i) = j. To visualize 
this definition differently, let ax and ay be the permutations of order n such that < 

Xa,{ 2 ) < • • • Xa,,(n] and < • • • Vayin), respectively (since the marginals of p are 
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uniform, ties do not occur with probability 1). Then the above definiton is equivalent to 
setting Tin = (Ty^ o ax- It is easy to see that if ^ has density p, then for any permutation 
Tin £ Sn one has 

„ n 

lPn« = TTn) = n! / p((i) j rfUjdUj. (1.7) 

By [11, Lemma 4.2] it follows that 'Kn converges weakly to p in probability. 

Our next corollary derives limiting distributions for p random permutations, when the 
measure p has a continuous density function with respect to Lebesgue measure. 

Corollary 1.7. Suppose iTn is a pp random permutation in Sn for some p G C. In this case 
the following conclusions hold: 

(a) If {crn}n>i is a sequence of non random permutations with an G Sn which converges to 
p, then Nn{TTn,o'n) Converges to Poi(p[p]) in distribution and in moments. 

(b) |c'„(l),--- ,Cnil)} converges to |Poi(cp(l)), • • • ,Po7(cp(/))| in distribution and in 
moments, where {Poi{cp{i))}\^^ are mutually independent. 

Even though the weak convergence of the random permutation sequence is the main ingredient 
in all the above results, the equi-continuity in both co-ordinates is not just a technical requirement. 
The following example shows that the conclusions of Theorems 1.1 and 1.3 might not hold if the 
equi-continuity condition fails. 

Proposition 1.8. Let be a probability distributon on Sn with the p.m.f. 

IRn,e(7rn) = 

where Cn is the identity permutation, and Nn{T^men) is the number of fixed points in Hn- Then for 
every 0 0 the following conclusions hold: 

(a) The random variable A^„(7r„,e„) converges to a Poisson random variable with mean e^ in 
distribution and in moments. 

(b) TTn converges weakly to u, the uniform distribution on [0,1]^ which is free of 6. 

(c) 

]Rn,e(7In(l) = l,7In(2) = 2) _ ^26 , 

Kn,e(7rn(l) = 2,7r„(2) = 1) 

Remark 1.3. Thus even though the sequence of random permutations under converge to 
Lebesgue measure (which is free of 9 and has a continuous density), the number of fixed points 
has a limiting Poisson distribution which depends on 9. This is the case as equi-continuity in both 
coordinates does not hold here, as demonstrated by part (c) of the proposition. 

1.2. Scope of future research. For the Mallows model with Kendall’s Tau, the results of this 
paper only apply for the case n(l — q{n)) = 0(1). If n(l — q{n)) —>■ oo, one should expect the 
number of fixed points to go to -|-oo, and computing the weak limits/limiting distribution after 
centering/scaling in this case remain open. In another vein, one might expect that convergence 
in the sense of permutations along with “mild” regularity conditions imply the weak convergence 
of LIS, as worked out for the Mallows model with Kendall’s Tau in [14]. Finally, computing the 
limiting density for the model defined in (1.6) might help give a more explicit description for the 
parameters of the limiting distributions of Corollary 1.6, as well as give non trivial copulas (bivariate 
distributions with uniform marginals) which constitute a subject area of its own in Finance. 
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1.3. Outline of the paper. Section 2 gives the proof of Theorem 1.1, Corollary 1.2, and Theorems 
1.3 and 1.4. Section 3 concludes the paper by proving Corollaries 1.5-1.7, and Proposition 1.8. 

2. Proofs of main results 

2.1. Proof of Theorem 1.1 and Corollary 1.2. 

Proof of Theorem 1.1. For A: G N setting 

i(7rn(p) = q) _ ^ 

ri(7r„(r) = s) 

condition (1.4) can be stated as 


en{k) := sup 

p,q,r,secS(n,Z) :||p-r||cx> <n/k 


lim lim en{k) = 0. 

k^oo n^oo 


( 2 . 1 ) 


Fix /c G N and partition (0,1] as with la := f ■ Setting 

i i 

a=l a=l 

note that ^p G A^, G Bif. Now for any r, s G S{n, 1) such that ^r G G Bif we have 

'n(7rn(p) = q) 


-1 


< en{k), 


IF’n(7rn(r) = s) 

which on summing over r G A^ , s G B^ and noting that the number of terms summed is at least 
(n — l)^^fc“^^ gives 

k'ii 

i(7rn(p) = q) <(1 + en(A:)) 7 y-^ Pn(7r„(r) = s) 


(n- 1) 


— (1 + ^nik)) 


k‘^^n‘ 


{n — 1)2* 

l * 


reAfe seRfc 




tv 

— (1 + ^n{k)) / _ 1 \2l 11 ^'^rXk\kpa/n\ ^ \kqa/n\\i 

a=l 


where denotes the I fold product measure of Using the fact that p is the density for pp 
this readily gives 


n¥„(7rn(p) = q) 
p,qe-S(n,z) nL=iP(^>lG) 


sup 


^(1 + (-nik)) X 


n 


21 


(n- 1) 


21 


xE sup 

ije[fc]' ni=l ^ 


sup 


\{\=lP{Xa-,Va) 


( 2 . 2 ) 

(2.3) 

(2.4) 


x,y,z,we[0,l]L||x-z||oo<l/Zc,||y-w||oo<l/Zc na=l P(^<o '^’a) 

The term in the r.h.s. of (2.2) converges to 1 on letting n —)• oo followed by A —>■ oo, using (2.1). 
Since U 7 r„ converges to pp, by [11, Theorem 5.2] we have 


max max 
ae[Z] ia&[k] 


^TTn [kia ^ ] Pp [Ip X Ij^ ] 


0 , 
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which along with the observation that x Ij^] is uniformly bounded away from 0 gives 

na=l ^ ^ja] P. 


max max 


1 


ae[Z] iam Hp[Ii, X IjJ 

as n —oo, for k fixed. An application of Dominated Convergence theorem implies that the term in 
(2.3) converges to 1 as well. Finally (2.4) is free of n, and converges to 1 as A; —)• oo by continuity 
of p. Combining this gives 

n¥n(vrn(p) = q) 


lim sup lim sup sup 

k^oo n^oo p,q,re5(n,0:||p-r||oo<n/fc Wa=lP{^^^) 


< 1, 


thus giving the upper bound in (1.5). A similar proof gives the lower bound, thus completing the 
proof of the theorem. □ 

We now introduce some auxiliary variables, to be used in the proofs of Corollary 1.2, and Theo¬ 
rems 1.3 and 1.4. 


Definition 2.1. For every n > 1 let {Z„(l), • • • , Z„(n)} be mutually independent random variables 
supported on [n] such that the marginal laws are given by 


for some p € C. Also set 


and for I > 1 set 


Qn{Zn{p) = q) 


p{pln,qln) 

YTs=iP{p/n,q/n) 


Mn{(7n) ■■= ^ HZnip) = <T„(p)}, 

p=l 

Dn{l):= ^ l{Z„(p)=r(p)}. 

pGU{n,l) 


Proof of Corollary 1.2. With as constructed in definition 2.1 we have 


|lPn,(7rn(Pn) = Q) “ Qn(^n(Pn) = q) 

(lGS{n,l) 


< max 
p,qeiS(n,Z) 


^nj'^njPn) = q) 
Qn{Z{pn) = q) 



Qn{ZniPn) = q) 

qeiS(n,Z) 


= max 
p,qeiS(n,Z) 


IP’n(7rn(Pn) = q) 
Qn{Zn{Pn) = q) 


which goes to 0 by (1.5). This implies that the laws of 7r„(p„) and Z„(p„) are close in total 
variation. Since the desired conclusion can be verified easily for Z„(p„), the proof is complete. □ 


2.2. Proofs of Theorem 1.3 and 1.4. We will use Stein’s method based on dependency graphs 
to prove Poisson limit theorems, as explained below: 

Let {Xa}a£i be a finite set of Bernoulli random variables. A dependency graph for {Xa}a^i is 
a graph with node set I and edge set E, such that if Ii, I 2 are disjoint subsets of I with no edges 
connecting them, then {Xa}a&h {^/ 3 }/ 3 e /2 independent. Let Na be the neighborhood 
of vertex a , i.e. A^(a) := {/3 G I : { oi ,( 5 ) G E } U {a}. Then one has the following Poisson 
approximation result, first proved in [1]. 
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Theorem 2.1. [7, Theorem 15] Let {Xa}a£i be a finite set of Bernoulli random variables with 
dependeney graph {I,E). Then setting A := YlaeiP<^’ ^ have 

||£(H/)-£(Poi(A))||Ty<X: E + ^ WiXffiEXp. 

a£l l3GN{a)/{a} a£l f3GN{a) 

The following lemma uses Theorem 2.1 to prove two Poisson limits which will be used in the 
proofs of Theorems 1.3 and 1.4. 


Lemma 2.1. Let Mn{an) and Dn{l) be as in definition 2.1. 

(a) If On converges to p, £ A4 in the sense of permutation limits, then we have M„(cr„) -4 Pfj,[p], 
and 

lim KMn{an)^ = 'EPoi{p,[p])^, for all k £ N. 

n^oo 

(b) For any I £ N we have Dn{l) Pcp{i)j ^.nd 

lim EDn{lf = EPoi(cp(/))^ for all k£N. 

n^oo ^ 


Proof. Setting m := info<a;,j;<i p{x,y), M := supg< 3 , y<]^ p{x,y) we have 0 < m < M < oo. 

(a) Since the random variables Xp = l{Zn{p) = crn{p)} for p = 1, 2, • • • , n are mutually inde¬ 
pendent, the dependency graph of {Xi,X 2 , ■ ■ ■ ,Xn} is empty. It then follows by Theorem 
2.1 that 


\\C{Mn{an)) - C{Poi{Xn))\\TV < 


n 


p{pln,an{p)ln) 

Eg=i Pip/n, q/n) 


2 


< 


1 m2 

n 


where 



p{pln,an{p)ln) 

Eg=i Pip/n, q/n) 


n^oo 



p{x,y)dp = p[p], 


and so Mn{an) converges to Poi{p[p\) in distribution. To conclude convergence in moments 
it suffices to show that limsup„_,.oo EMn{(Tn)^ < oo for every k £ N. To see this, set 


5(n, 1) := {p £ S{n, 1) : Pi < p 2 <■■■ < pi} 


denote the set of all n tuples in increasing order, and note that 

k 

EMn{an)^ = ^ Qn{Zn{p) = CTnip)) < ^ Qn{Zn{p) = O-nCp))- 

pe[n]'' 1=1 p£S{n,l) 


Here the factor k^ in the r.h.s. above accounts for the fact that a specific term {Zn{p) = 
<7n(p)} with p £ S{n,l) can arise from at most terms in [n]^. Since |5(n, Z)| = (”), we 
can bound the r.h.s. above by 


k I 

E^' E n 


pSiS(n,Z) 


p{pa/n,an{pa)/n) ^ 
YTci,=iP^Paln,qaln) “ 



< oo. 


(b) The proof of part (b) is similar to the proof of part (a). For p G lA{n,l) setting Xp = 
l{Zn(p) = T(p)} note that Xp is independent of Xq whenever the indices p and q have no 
overlap. Thus the dependency graph of the random variables {Xp, p £ lAn,i} has maximum 
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degree at most Also for any p, q which overlap we have EXpAq = 0 unless p = q. 

Thus an application of Theorem 2.1 gives 

f 7l\ f T) — l\ 1 

||£(O„(0) - £(P.JII S (,) (i - 1)! X (, _ Ji! X - X 

with 

, ^ 1 PiPiln,P2ln) p{Piln,piln) n^oo 

” " ~l E;=i p{pi/n. qi/n) ^ ‘ ‘ ‘ ^ ^^=1 P{Pl/n, qi/n) ^ 

and so Dn{l) converges to Poi{c{l)) in distribution. Convergence in moments follows by a 
similar calculation as before. □ 

Proof of Theorem 1.3. Let {Z„(l), • • • ,Z„(n)} and Mn{an) be as defined in 2.1. Then using part 
(a) of Lemma 2.1 and the fact that the Poisson distribution is characterized by its moments, it 
suffices to show that for every A: G N we have 

lim |EAn(7rn,cr„)*'- EMn(c7n)*'| = 0. 

n^oo 

To this effect setting Zn(p) = (Zn(pi), ■ ■ ■ , Zn(pk)) for p G [n]^ we have 

|EAn(7r„,C7n)*' - EM„((Tn)^| < ^ ||Pn(vrn(p) = O'n(p)) -Qn(^-^n(p) = 0-„(p)^ | 

pS[n]'' 

First note that the events {vr„(p) = cr„,(p)} and {Zn(p) = (Tn(p)} have positive probability for all 
p e [n]^, and so for any p G [n]^ setting L = L{p) denote the number of distinct indices gives the 
bound 


|Pn(7r„(p) = an(p)^ - Qn(^Zn{p) = O'n(p)) | 


< max 
p,qeiS(n,L) 


(7rn(p) = q) 


i{Znip) = q) 


- 1 


Zn{p) = O'n(p))- 


Since L(p) < k, taking a maximum over L and summing over p G [n]^ gives the bound 


|EA„(7r„,a„)'^ 




< max max 

[ l^[k] p,qG(5(n,/) 


IP’n(7rn(p) = q) 
Qn{Zn{p) = q) 



EMnian)’^. 


By (1.5) we have 


Since Lemma 2.1 implies 


max max 
lG[k] p,qecS(n,Z) 


IPn(7rn(p) = q) 
Qn{Zn{p) = q) 


limsupEMji(cr„)^ = EPoi(/i[p])^ < oo, 

n^oo 


the r.h.s. of (2.5) converges to 0 as n —)• oo, thus completing the proof of the theorem. 


(2.5) 


□ 


Proof of Theorem l.f. Let {Z„(l), • • • , Z„(n)} and {Il„(o), 1 < a < /} be as defined in 2.1. Then 
by part (b) of Lemma 2.1, for any finite collection of non negative integers ki,k 2 , - ■ ■ ,ki we have 

i i 

lim TT ET>„(a)^“ = TT EPoi(cp(a))’^P 

a=l a=l 
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Thus to complete the proof it suffices to show the following: 

i i 


lim 

n^oo 


lim 

n^oo 


'n[a) - 

a=l a=l 

I I 




a=l 


a=l 


= 0 , 


= 0 . 


( 2 . 6 ) 

(2.7) 


For showing (2.6) we have 

i 


a=l 


a=l 


|E [] Z7„(a)'^“ -llEDn{a)’^‘^\ |Qn(nti n^;^,{z„(p(a, 6,)) = T(p(a,6a 

r 

i 

n‘Q^(^fer=i [Zn{p{a,ba)) = T{p{a,b. 


a=l 


( 2 . 8 ) 


where 

r := |p(a, ba) £ U{n, a), ba = 1,2, ■■ ■ , ka, a = 1,2, ■■ ■ , /|. 

Proceeding to analyze a generic term in the r.h.s. of (2.8), fix 

p{a,ba) £U{n,a), 1 < ba < ka,l < a < 1. 

Let La = La{p{a,ba),l < ba < ka} denote the set of distinct indices in the set {p(a,6a),l< 
ba < ka}- First note that if the sets La do not overlap across a, both terms in the r.h.s. of (2.8) 
are the same, and so gets canceled. As an example, this happens for the choice 


I = 3,k^ = 0,k2 = l,k5 = 2, p(2,l) = (1,2), p(3,l) = (3,4,5), p(3, 2) = (3,4,5). 

In this case Li = cf), L 2 = {1,2} and L 3 = {3,4, 5} do not overlap, and so the corresponding terms 
in the r.h.s. of ( 2 . 8 ) get cancelled. 

If the sets La do overlap across a, then the first term in the r.h.s. of (2.8) is 0. In this case 
setting L := X)a=i total contribution of the second term in the r.h.s. of (2.8) is bounded 

by (^)^. Since there is a repetition among the indices, the number of distinct indices L{D) in 

the set {p(a, 6 a), 1 < 6 a < fea, 1 < a < f} is strictly less than L. As an example, this happens for 
the choice 


I = 3,ki = 0,k2 = l,k3 = 2, p(2,l) = (1,2), p(3,l) = (3,5,4), p(3, 2) = (1, 6 , 7). 


In this case Li = 0, L 2 = { 1 , 2 },L 3 = {1,3,4,5, 6 ,7}, and so the number of distinct indices 
L{D) = 7 which is less than L = IL 2 I + IL 3 I = 8 . Setting K := Yl\=ikai the total number of 
terms with exactly L{D) distinct indices is at most Summing over the possible ranges 

L[D) G [1, L — 1], L G [1,K] the total contribution of such terms is at most 


K L-l 

E E 

L=1L(D)=1 


n 

L{D) 


K\ 


M y 
mnJ 


= O 


thus proving ( 2 . 6 ). 
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Proceeding to prove (2.7) we again have 


i 

a=l 


I 

E H 

a=l 


- X] ( ^“=1 ^K=i'^n{v{a, ba)) = T{p{a, ba))j 
r 

~Qn ^ ba)) = T(p(a, ba))'j 


(2.9) 


Proceeding to bound the r.h.s. of (2.9), note that in this case if all the indices in the set 

are not distinct (i.e. L[D) ^ L), then both terms in the r.h.s. of (2.9) are 0. Even if L{D) = L, it 

is possible that both terms are 0, which happens for example for the choice 


I = 3,ki = 0,k2 = l,k3 = 2, p(2,1) 


(1,2), p(3,l) = (3,5,4), p(3,2) 


(3,4,5). 


In this case Li = {1, 2}, L 2 = {3,4, 5} and so L{D) = L = 5. However both the terms on the r.h.s. 
of (2.9) have 0 probability. If either of the terms have non zero probability, then a generic term 
on the r.h.s. of (2.7) is of the form |P„(7r„(p) = q) — QniZnip) = q)| for some p, q G S{n,l) with 
I £ [L], Noting that L < K, this can be bounded by 


max max 
1&[K] p,qg5(n,Z) 


En(7rn(p) = q) 
Qn{Zn{p) = q) 


nL=i n,": 


=1 




(p(a,6a)) = r(p(a,6a))) 


On summing over P using (2.9) gives 


IE TT C'n(a)^“ — E TT Dn(a)^°- \ < max max 

ielKW,ceSin,i) 

from which (2.7) follows on using (1.5) along with (2.6). 


i(7rTi(p) = q) 


i{Zn{p) = q) 


i|e Dn{c 


a=l 


□ 


3 . Proof of Corollaries 1.5-1.7 and Proposition 1.8 


Proof of Corollary 1.5. By [17, Theorem 1] it follows that 7r„ converges weakly in probability to 
the measure fip^ induced by the density defined in (1.2). Given Theorems 1.1, 1.3 and 1.4, for 
proving both parts (a) and (b) it suffices to verify the equi-continuity condition (1.4), which is 
equivalent to the following two conditions: 


lim lim sup 

5—>-0 n^oo p^q_rg5(n,Z):||p—r||cx><n(5 

lim lim sup 

(5—^-O n—>-oo q^r,seiS(n,Z):||q—s||oo<n(S 


En(7rn(p) = q) 
En(7rn(r) = q) 

En(Tn(r) = q) 

En(7rn(r) = s) 



= 0 , 

= 0 . 


(3.1) 

(3.2) 


Recall that (3.1) was already verified in [3, Corollary 6.3-|-Lemma 7.1]. By repeating the argument 
presented there, we prove both (3.1) and (3.2) here for completeness. To show (3.1), fix p, q, r such 
that Up — rjj^ < n6. Let H(p,q) denote the set of all permutations in Sn such that 7rn(p) = q, 
and H(r, q) be defined likewise. We will now define a bijection <I> = 4>[(p, q); (r, q)] from H(p, q) to 
H(r, q). For any 7r„ G H(p,q) set 

4'(7rn)(r) = q, 4>(7r„)(p) := 7r„(r), 4>(vr„)(i) = 7r„(i) otherwise. 


It is easy to see that <I> is indeed a bijection, and 


^n,q{n) ij^n) 




nlS 
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where we use the fact that the inversion status of a pair (i,j) in vr^ is the same as its inversion 
status in <h(7r„) unless i G '-J‘a=i\Pa.^ ^a] 3 ^ Q- Summing over 7r„ G f^(p, q) gives 


i(7rn(p) = q) 


< max 


{q{n),q{n) ^ 


nl5 


^^(TrnCr) = q) 

and since the bound in the r.h.s. above is free of p, q, r, taking a sup gives 


sup 


i(7rn(p) = q) 


< max 


{q{n),q{n] 


-1 


nlS 


p,q,rSiS(n,Z):||p—r||cx><fi<5 q) 

On letting n —>■ oo followed by 5 —)• 0 and noting that n(l — q{n)) —^/3 g (— oo, oo), we get 

i(7I'n(p) = q) 

(5—>-0 n—>-oo p,q,rG(S(n,/):||p— r||oo<r 


limsuplimsup sup '"i ^ < 1, 

,<ns ^nMr) = q) 


thus giving the upper bound in (3.1). By symmetry we have 

PnCvrnCp) = q) 


lim inf lim inf 


................ sup m, / / \ 

5^0 n^oc p,q,re-S(n,Z):||p-r|U<n5 ^nCTTnlr) = q) 


> 1, 


thus giving the lower bound, and hence proving (3.1). For proving (3.2) a similar argument works, 
except now we set up the bijection = <h„[(r,q); (r, s)] between Or,q to flr,s by setting 

$(7rn)(r) = s, $(7r„)(7r“^s) := q, $(7rn)(z) = 7rn(z) otherwise. 


The rest of the argument repeats itself, and we omit the details. 


□ 


Proof of Corollary 1.6. (a) It follows from [15, Theorem 1.4] that vr^ converges to a unique 

measure yiffi weakly in probability, which is the solution of the optimization problem 

where u is the uniform measure on the unit square, and T)(.||.) is the Kullback Leibler 
divergence. It was further shown there that yLffi has a density of the form gffi{x,y) = 
g^fi^,y)+ 0 'f, 9 {^)+bf,eiy)^ where af^Q{.) and bf^ 0 {.) are unique almost surely. To complete the 
proof of part (a), it suffices to show that the function gj^Q is continuous on the unit square, or 
equivalently that is continuous. To this effect, using the fact that has uniform 

marginals, we have 

g-a/,0(x) _ f ^Sf(x,y)+bf^g{y)^y^ 


which readily gives 




for almost all x G [0,1], and consequently is integrable. But then we have 

r-l 


,-a/,9(a:i) _ p-affi(x2) 


< sup 
y6[0,l] 


gSf{^i,y) _ Q6f{x2,y) 


hAy)dy, 


from which continuity of e follows from continuity of /(.,.). 
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(b),(c) As in the proof of Corollary 1.5 it suffices to verify the conditions (3.1) and (3.2). Using 
the same notations as in the proof of Corollary 1.5, we have 

Qn,f,e{'^n) _ fiPaln,qaln)-f {rain,qain) 


Qn,f,ei^{Trn)) 
and the exponent in the r.h.s. above is bounded by 

|6'| sup \f{xi,y)-f{x 2 ,y)\. 

xi,3:2,y&[0,l]:\xi—X2\<S 

Since this goes to 0 as (5 0, a similar proof as before verifies (3.1). The proof of (3.2) is 

similar, and again we omit the details. □ 

Proof of Corollary 1.7. Since a sequence of y,p random permutations converge to fip weakly in 
probability, it suffices to verify (3.1) and (3.2). 

To this effect, with (Xi, Yi), • • • , {Xn, Yn) first note that marginally both (Xi, • • • , Xn) 

and (Yi, • • • , Y„) are i.i.d. 17(0,1). Thus if {Ui, • • • , Un) and (Yi, • • • , Yi) are the order statistics of 
(Ai, • • • , Xn) and (Yi, • • • , Y„) respectively, for any <5 > 0 we have 


Ui-- 

n 


> -f J + Pn —h <5^ 

by Hoeffding’s inequality. Also using (1.7), for any p,q G S{n,l) we have 


< z < 2e 


—S^r 


(3.3) 


li-^nip) = q) =nl 


7r„efl(p,q) ,Ua,vi<V2<-<v„ 


'[{f[ui,v^^{i)]duidvi 


which, for r € S{n, 1) gives 


i(7rn(p) = q) T-r p{xa,ya) /MY 

^(7^n(p) = r) “ x,y,ze[0,l]'f=i ' 


Noting that |5(n, 1)| < n\ summing over r this implies 

/ m C 


(j^) S ='^) = (;s) 


reiS(n,Z) 

Finally, for any p,q,r G S{n,l) such that ||p — r||oo < n5, setting An := {maxjgj^] 
we have 

Pn(vrn(p) = q) ^Pn(^n) +Pn(Tn(p) = q, A„) 


Ui-f 


(3.4) 


<5} 


i(7r„(r) = q) 


Pn(Tn(r) = q, A„) 

Y n- X/ fui<U2<-"<u„,vi<V2< 
7r„en(p,q) 


••• ,<v„,A„ ni=l '^7r„(i)) 


n: 


E L 


7r„er2(r,q) 


Ul<U2<---<Un,Vi<V2<-" ,<V, 


,A„ ni=l p{f^iiX-Kn{i^ 


< 


\{\=lp(upa,V^ 


max 


sup 


■'Ha 


+ 


p.q,re5(n.Z):||p-r||^<n5^_^g[0_l]„^|„,_l|<^ IlLl P 
Pn(vrn(r) = q) - Pn(A^)' 


(3.5) 


(3.6) 
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Since \up^ — UrJ <25+ < 3(5, the expression in (3.5) can be bounded by 

W\=l pi.Xa,ya) 
x,y,ze[0,l]*:||x-z||oo<35 na=l Pi^a^ Va) 

which is free of n, and goes to 0 as (5 —)■ 0 by continuity of p. Also, using (3.3) and (3.4) it follows 
that the expression in (3.6) is bounded above by 

2g-n52 

~/ W ’ 

I ^ \ _ 

I nm J 

which converges to 0 as n —>■ oo, for every 5 hxed. Thus, taking a maximum over p, q, r G S{n, 1) 
such that Up — r|loo < nS we have 

1- IF’n(vrn(p) = q) / 1 

hmsuphmsup max ——— -- < 1, 

5-5.0 n^oo p,q,re*S(n,i),||p-r||oo<n(5 r„(7r„(r) = q) 

thus giving the upper bound in (3.1). Similar arguments give the lower bound in (3.1), as well as 
(3.2), thus completing the proof of the corollary. □ 


Proof of Proposition 1.8. With = Mn,o denoting the uniform measure on Sn and denoting 
the number of derangements of n, we have 


_ jgp g6Af„(7r„,e„) 

n! ”■ 


E' 

/c=0 


Ok 




n\ 


exp{e® - 1}, 


(3.7) 


where we use the fact that Dn/n\ converges to e 

(a) For any A > 0 we have 

^ exp{e®(e^ - 1)}, 

and so Nn{TTn+n) Converges to Poi{e^) in distribution and in moments. 

(b) With i7(.||.) denoting the Kullback-Leibler divergence we have 

/g.z„(0). 

D(M„_o||Kn, 0 ) = log [ - 6'IEp„A^(vrn,en) e® - 1 - 0, 

and so by [4, Prop 5.1] we have that the two probability distributions g and Mn,o = IPn 
are mutually contiguous. Since 7r„ converges weakly to u under = Mri, 0 ) by contiguity 
the same happens for 

(c) Let An := {vr^ G Sn : vr„(l) = l,7r„(2) = 2}, and Bn := {vr„ G Sn ■ vr„(l) = 2,7r„(2) = 1}. 
Define a bijection oj from An to Bn by setting a;(7r„)(i) = i for 3 < i < n, and note that 

l^n,6('^n) _ 26 

Kn,6»(w(7rn)) 

and so summing over vr^ G An gives 

l^n,6(7I'n(l) = 1) 71n(2) = 2) _ 20 / 

M„,0(7r„(l) = 2,7r„(2) = l) ^ ’ 


thus proving part (c). 


□ 










FIXED POINTS AND CYCLE STRUCTURE OF RANDOM PERMUTATIONS 


17 


4. Acknowledgements 

The Poisson distribution for the number of fixed points in the Mallows model with Kendall’s 
Tau was conjectured by Susan Holmes based on empirical evidence. This paper also benefited from 
helpful discussions with Shannon Starr. Suggestions from an anonymous referee greatly improved 
the presentation of the paper. 


References 

[1] Arratia, R., Goldstein, L., and Gordon, L. (1990). Poisson approximation and the Ghen-Stein method. 
Statist. Sci. 5, 4, 403-434. With comments and a rejoinder by the authors. MR 1092983 

[2] Basu, R. and Bhatnagar, N. (2016). Limit Theorems for Longest Monotone Subsequences in Random Mallows 
Permutations. Available at http://arxiv.org/pdf/1601.02003. 

[3] Bhattachara, B. and Mukherjee, S. (2015). Degree sequence of random permutation graphs. Ann. Appl. 
Probab., to appear. 

[4] Bhattachara, B. and Mukherjee, S. (2015). Inference in Ising models. Available at http://arxiv.org/abs/ 
1507.07055. 

[5] Bhatnagar, N. and Peled, R. (2015). Lengths of monotone snbsequences in a Mallows permutation. Probab. 
Theory Related Fields 161, 3-4, 719-780. MR 3334280 

[6] Borodin, A., Diaconis, P., and Fulman, J. (2010). On adding a list of numbers (and other one-dependent 
determinantal processes). Bull. Amer. Math. Soc. (N.S.) 47, 4, 639-670. MR 2721041 

[7] Ghatterjee, S., Diaconis, P., and Meckes, E. (2005). Exchangeable pairs and Poisson approximation. 
Probab. Surv. 2, 64-106. MR 2121796 

[8] Diaconis, P. (1988). Group representations in probability and statistics. Institnte of Mathematical Statistics 
Lecture Notes—Monograph Series, 11. Institute of Mathematical Statistics, Hayward, CA. MR 0964069 

[9] Diaconis, P. and Ram, A. (2000) Analysis of Systematic Scan Metropolis Algorithms Using Iwahori-Hecke 
Algebra Techniques. Michigan Math. J. 48, 1, 157-190. MR 1786485 

[10] Gladkich, a. and Peled, R. (2016) On the cycle structure of Mallows permutations. Available at http: 
//arxiv.org/pdf/1601.06991. 

[11] Hoppen, C., Kohayakawa, Y., Moreira, C. G., Rath, B., and Menezes Sampaio, R. (2013). Limits of 
permutation sequences. J. Combm. Theory Ser. B 103, 1, 93-113. MR 2995721 

[12] Kenyon, R., Kral, D., Radin, C., and Winkler, P.(2015). A variational principle for permutations. Avail¬ 
able at http://arxiv.org/pdf/1506.02340. 

[13] Mallows, G. L. (1957). Non-null ranking models. 1. Biometrika 44-: 114-130. MR 0087267 

[14] Mueller, G. and Starr, S. (2013). The length of the longest increasing subsequence of a random Mallows 
permutation. J. Theoret. Probab. 26, 2, 514-540. MR 3055815 

[15] Mukherjee, S. (2016). Estimation in exponential families on permutations. Ann. Statist. 44, 2, 853-875. 
MR 3476619 

[16] Nelsen, R. B. (2006). An introduction to copulas, Second ed. Springer Series in Statistics. Springer, New York. 
MR 2197664 

[17] Starr, S. (2009). Thermodynamic limit for the Mallows model on Sn. J. Math. Phys. 50, 9, 095208, 15. 
MR 2566888 

[18] Walters, M. and Starr, S. (2015). A note on mixed matrix moments for the complex Ginibre ensemble. J. 
Math. Phys. 56, 1, 013301, 20. MR 3390837 

[19] Trashorras, j. (2008). Large deviations for symmetrised empirical measures. J. Theoret. Probab. 21, 2, 397- 
412. MR 2391251 

Department of Statistics, Golumbia University, New York, USA, sm3949@columbia.edu 


